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Abstract: Dijkgraaf and Vafa (DV) have conjectured that the exact superpotential 
for a large class of M = 1 SUSY gauge theories can be extracted from the planar 
limit of a certain holomorphic matrix integral. We test their proposal against existing 
knowledge for a family of deformations of M = 4 SUSY Yang-Mills theory involving 
an arbitrary polynomial superpotential for one of the three adjoint chiral superfields. 
Specifically, we compare the DV prediction for these models with earlier results based 
on the connection between SUSY gauge theories and integrable systems. We find 
complete agreement between the two approaches. In particular we show how the DV 
proposal allows the extraction of the exact eigenvalues of the adjoint scalar in the 
confining vacuum and hence computes all related condensates of the finite- N gauge 
theory. We extend these results to include Leigh-Strassler deformations of the M — 4 
theory. 



1. Introduction 

Theories with Af = 1 supersymmetry in four dimensions are of considerable interest. 
Apart from their potential phenomenological relevance, they provide a field theoretic 
laboratory where physical phenomena such as confinement and dynamical symmetry 
breaking can be studied in a controlled setting. This is possible largely because these 
theories contain a special class of observables whose dependence on the parameters 
of the theory is constrained to be holomorphic. These quantities, which include the 
vacuum expectation values of chiral operators and the tensions of BPS saturated 
domain walls, are governed by a holomorphic superpotential which can sometimes 
be determined exactly. Although exact results for many different models have been 
obtained, a general recipe for computing the superpotential for an arbitrary Af = 1 
model has remained elusive. In recent work [1-4], Dijkgraaf and Vafa (DV) have 
proposed a solution to this problem which applies to a large class of Af = 1 theories. 
Specifically, they claim that the superpotentials for these theories can be extracted 
from the planar limit of a certain holomorphic matrix integral. In this paper, we 
will perform a detailed test of their proposal against known results for a family of 
Af = 1 preserving deformations of Af = 4 SUSY Yang-Mills theory with gauge group 
SU(N) (or U(N)). For this family of theories there is an alternative method for 
computing the exact superpotential, given in [5-8]. In fact, the superpotential can 
be deduced from the value of the Lax matrix of a certain classical integrable system 
in its equilibrium position. Remarkably, we will find that these two seemingly very 
different recipes give exactly the same result. The comparison will also shed light on 
both approaches. In the remainder of this introductory section, we will review both 
the DV proposal and the results of [8] based on the connection to integrable systems 
and present our main conclusions. 

We begin by considering Af = 4 SUSY Yang-Mills theory with gauge group G = 
SU(N) and bare complex coupling r = Ani/ gy M + 9 /2-n . In the standard language of 
Af = 1 superfields, the theory contains three adjoint- valued chiral multiplets, which 
we will label $ + , $~ and $. We will consider the deformation of the Af = 4 theory 
obtained by adding an Af = 2 preserving mass term for $ ± , together with a general 
polynomial superpotential for <£> which preserves Af — 1 supersymmetry. The classical 
superpotential then reads, 




(1.1) 



Setting g p = for p > 3 corresponds to the Af = 1* deformation but we will not 
limit ourselves to this case. 
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This theory already has a rich vacuum structure at the classical level. In partic- 
ular, there is a large set of classical vacua in which the adjoint scalar fields acquire 
vacuum expectation values (VEVs) spontaneously breaking some or all of the gauge 
symmetry. In this paper we will focus on SUSY vacua in which the full SU(N) gauge 
symmetry is unbroken. At the classical level this corresponds to setting the adjoint 
scalar VEVs to zero 1 . The qualitative behaviour of the corresponding quantum the- 
ory in these vacua is well known: at energy scales much less than the mass scale set by 
the deformations, the chiral multiplets decouple leaving the theory of a single SU(N) 
vector multiplet which is in the same universality class as M = 1 SUSY Yang-Mills. 
The theory is therefore in a confining phase with a mass gap and gluino condensation. 
There are N supersymmetric vacua denoted \k) with k = 0, 1, 2, . . . , N — 1, each with 
a non-zero expectation value for the gluino bilinear which is the lowest component of 
the superfield S = Tr (W a W a )/327i 2 . In fact, the action of S-duality, inherited from 
the underlying M = 4 theory, dictates that the vacuum states \k) and \k + 1) are 
simply related by a 2ir shift in the vacuum angle: r — ► r+ 1. Hence it will suffice to 
consider the k = vacuum, |0), in the following. Our aim will be to determine the 
exact effective superpotential in this vacuum, (OlWcfjlO). The vacuum expectation 
values of the operators u p = Tr $ p can then be obtained by differentiating with 
respect to the parameters g p appearing in (|I~T|). 

1.1 The integrable system approach 

An exact formula for the effective superpotential can be obtained from a slight 
extension of the results presented in [8]. This approach relies on the connection 
between SUSY gauge theories and integrable systems first discovered in [9-11]. In 
the present case the relevant integrable system is the elliptic Calogero-Moser (ECM) 
system. This is a system of N non-relativistic particles with complex positions and 
momenta, Xj and Pj, j = 1, 2, . . . , N, which live on a torus E T of complex structure 
r and its tangent space respectively. Time evolution is governed by the quadratic 
Hamiltonian, 

N p2 

3=1 j>k 

where p is the Weierstrass function (see Appendix A.l) which is elliptic with respect 
to the torus E T . The system has a Lax formulation involving two N x N matrices 
and Mjk, in terms of which Hamilton's equations can be rewritten as iL = [M,L]. 
This immediately implies the existence of N conserved quantities H p = Tr L p for 

general, the adjoint scalar fields can and will acquire VEVs in the quantum theory without 
breaking the gauge symmetry. 
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p = 1, . . . , N. These quantities can also be thought of as a set of Poisson-commuting 
Hamiltonians each generating a distinct time-evolution of the system. 

We will now state the relation between the deformed M = 4 theory and the 
ECM theory (for more details see [8]). Firstly, for each value of p, the conserved 
quantity H p is identified with the expectation value of the operator u p = Tr $ p 
introduced above 2 . As explained in [6], the supersymmetric vacua of the theory then 
correspond to the equilibrium configurations of the ECM particles. In particular, 
for the theory with an arbitrary polynomial deformation, as in Eq. (|l.l| ) we must 
consider equilibrium configurations with respect to the time-evolution generated by 
the Hamiltonian H = J2p=2 9p^p- 

In the general case, finding these configurations is not straightforward. However, 
the massive vacua of the theory have a very special property: they correspond to 
static configurations which simultaneously stationarize each of the Hamiltonians 3 
Hp, for p = 2, . . . , N. In particular, the vacuum |0) introduced above corresponds to 
a configuration where the particles are evenly spaced around one cycle of the torus: 
Xj = (2irir/N)j - mr(N - l)/N; Pj = for j = 1, . . . , N. As each Hamiltonian is 
stationary, the same configuration corresponds to a SUSY vacuum for all values of the 
deformation parameters. The eigenvalues of the Lax matrix L a b in this equilibrium 
configuration were calculated explicitly in [8]. The result is 

~ 1 9' 3 (n/2 - nj/N\r/N) 

3 ~ 2 6 3 (7r/2-nj/N\T/N) 

(1.3) 

where j = 1,2, ...,N. The vacuum expectation values of the operators, u p , can 
then be determined via the relation (0|w p |0) = Ylf=i ^j- As ^ ne resulting formulae 
are completely independent of the deformation parameters g p , we find that that the 
vacuum value of the effective superpotential is simply given by, 

N N N 

(0\W eS \0) = 9 P (0K\0) = E E ^ (1-4) 

p=2 p=2 j=l 

2 To restrict to the SU(N) theory where u\ vanishes we consider the ECM system in its center- 
of-momentum frame where H\ = J2j=i Pj = 0- 

3 In the field theory this is related to the fact that massive vacua preserving N — 1 supersymmetry 
correspond to the maximal degenerations of the Riemann surface which governs the Coulomb branch 
of the Jsf = 2 theory obtained by setting all the deformation parameters to zero. However, it can 
also be demonstrated directly using the Lax representation for the equations of motion of the ECM 
system. This is shown in Appendix B. 
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1.2 The Dijkgraaf-Vafa Approach 



Dijkgraaf and Vafa have proposed an exact recipe for calculating the exact effective 
superpotential W e s as a function of the gluino condensate (the expectation value 
of the superfield S). Specifically, they relate the superpotential to the following 
holomorphic matrix integral: 

Z = J [c/$p$+p$-] exp f-yWi^ (1.5) 

where g s is a coupling constant. The integral is to be understood as a multi- variable 
contour integral. As for any adjoint SU(N) theory the free energy T has a 't Hooft 
expansion where the Feynman diagrams are organised according to their genus g > 0. 

where A = g s N is the matrix model 't Hooft coupling. DV predict that the exact 
superpotential of the above Af = 1 perturbation of the M = 4 theory only depends 
on the leading term in this expansion: J-'o(X). This term may be evaluated by 
taking the 't Hooft limit of the matrix model, N — > oo, g s — > with A held fixed, 
which isolates the planar Feynman diagrams. Fortunately the matrix model can 
be solved exactly in this limit using standard saddle-point technology. In the case 
of a quadratic deformation, the solution has been already been given by Kazakov, 
Kostov and Nekrasov [13]. Generalization of their solution to an arbitrary polynomial 
deformation is straightforward and is given below. 



The proposed dictionary between matrix model and field theory quantities identi- 
fies the 't Hooft coupling A = g s N with the expectation value of the gluino bilinear 
superfield S. From now on we will use the symbol S to denote both the matrix model 
't Hooft coupling and the gluino condensate of the field theory and assume that the 
interpretation is clear from the context. The resulting expression for the field theory 
superpotential is, 

W eS (S) = N-£ - 2mrS (1.7) 



The supersymmetric vacua may be found be extremising W e R as a function of S. 
The expectation value of W e g in each vacuum is simply its value at the corresponding 
the extremum. The main aim of this paper is to compare the resulting formulae with 
the exact field theory predictions Eqs. ( |1.3| ) and (|1.4|). We find that the DV prescrip- 
tion allows us to extract the field theory eigenvalues Xj and we find precise agreement 
between the two approaches, for arbitrary values of the deformation parameters g p , 
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up to the known ambiguities in the definition of the operators u, p [5]. Agreement in 
the case of a quadratic superpotential has previously been checked in [4]. 

Several features of our result are worthy of further comment. A remarkable 
consequence of Dijkgraaf and Vafa's proposal is that the holomorphic quantities of 
the deformed J\f — 4 theory are large- N exact. Although this property is far from 
obvious in the integrable systems approach it nevertheless seems to be true that 1/N 
corrections to all physical quantities vanish up to the usual operator mixing ambigu- 
ities! On another hand there are some features of the field theory results which are 
not apparent at the outset in the DV prescription, for example the linear dependence 
of the exact superpotential (|1.4j) on the deformation parameters g p . Nevertheless, we 
are able to prove that for the class of deformations considered the DV prescription 
automatically implies this linearity in the confining vacuum. 

We also extend our study to include a different class of deformations of the 
M = 4 theory, namely the Leigh-Strassler deformations. We use the DV proposal 
to compute the exact effective superpotential for these deformations and extract the 
eigenvalues Xj in the confining vacuum of the theory. Remarkably, we find once again 
that these are independent of the deformation parameters g p . 

2. Superpotentials and eigenvalue distributions from matrix 
models 

According to the Dijkgraaf-Vafa proposal [1-3] superpotentials for a large class of 
M = 1 supersymmetric field theories at finite N are computed by corresponding 
large-A" matrix models. Assuming the validity of this conjecture for arbitrary de- 
formations of the M = 4 theory, we will demonstrate that for such deformations, 
in the confining vacuum of the resulting theory the matrix model in fact yields the 
correct eigenvalue distribution for the adjoint scalar field in the theory. This tests 
the Dijkgraaf-Vafa proposal for a large class of deformations of the M = 4 the- 
ory. Furthermore, it provides a general prescription for extracting QFT eigenvalue 
distributions from the matrix model and also demonstrates that the matrix model 
computes all condensates in the holomorphic sector of the theory. 

2.1 Deformations of J\f = 4 theory 

In the M = 1 language, the M = 4 U(N) theory with coupling constant r = 
Airi/gyM + #/27r has three adjoint chiral superfields $ + ,$~ and $. We consider a 



5 



general class of deformations of the M = 4 theory specified by a tree level superpo- 
tential 

W = ^Tr [«$[$+, + +V(<$>)] (2.1) 

9ym 

where V($>) is a general polynomial 

vm^gp**- (2-2) 

p 

Note that for the sake of simplicity we have set any and all mass scales to unity. (For 
V(<&) = $ 2 we obtain the J\f = 1* theory.) 

As a direct consequence of the proposal of Dijkgraaf and Vafa [1-3], in a given 
vacuum the exact effective superpotential for the above class of deformations of the 
M = 4 theory is computed by the planar diagram expansion, i.e. a large- N limit, of 
the three-matrix model partition function expanded around that vacuum 

Z= /[d$ + p$-p$]exp-— Tr [i$ [$+,$"] + $ + $~ + 1/($)] (2.3) 
J 9s 

with $+ = ($-) t . 

The confining vacua of the field theory correspond to the classical solution $ = 
$ ± = which preserves the full U(N) gauge symmetry. In order to implement the 
Dijkgraaf- Vafa proposal for extracting the field theory superpotential in the confining 
vacua we therefore need to solve the matrix model Eq. ( |2.3j ) around the trivial vacuum 
$ = $ ± = 0. As noted for example in [13], the key fact that permits the solution of 
these matrix models is that one may integrate out $ ± exactly to obtain a one-matrix 
integral to be solved in the large- N limit 

Z= / d$ ±— -. (2.4) 

J 1 J det(Adj 4> + ^) v ; 

2.2 Solution of the large-A^ matrix model 

The above one-matrix model actually becomes tractable in the large-iV limit by 
going to an eigenvalue basis and performing a large- saddle-point approximation 
to the integral. The details of this procedure have been extensively discussed in the 
literature. We refer the reader to [13], and references therein, for details. 

The main feature of the solution is that the eigenvalues Xj of $ interact via a 
repulsive effective potential and form a continuum in the large- N limit and condense 
into a cut along the real axis. The actual extent of the cut and the density of 
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eigenvalues p(A) along the cut is self-consistently determined by the saddle point 
equation in terms of the parameters of the deformation V(<&) and the matrix model 
't Hooft coupling S = g s N. The saddle point equation in the large- N limit is most 
conveniently written in terms of the resolvent function 

oo(z)= [ + ^-dX AG [-a, a]; / + p(X)dX = 1. (2.5) 
J -a z ~ A J - a 

The resolvent u(z) is an analytic function on the complex z-plane whose only sin- 
gularity is a branch cut with the discontinuity across the cut given precisely by the 
spectral density 

uj(X + ie) -w(A-ie) = -2mp(X); Xe[-a,a]. (2.6) 



The saddle point equation essentially encodes the condition that the force on a 
test eigenvalue due to the remaining large- N equilibrium distribution of eigenvalues 
vanishes along the cut where the eigenvalues condense. The matrix model effective 
action yields this force on a test eigenvalue at position z in the complex plane as 

F(z) = ^[G(z + i)-G(z-i)] (2.7) 

where 

G(z) = U(z) + iS(uj(z + - uo{z - |)). (2.8) 
and U(z) is a polynomial in z such that 

V'{z) = \[U{z + ±)-U(z-i)]. (2.9) 



The analytic properties of the function G(z) are central to the solution of the 
model and its subsequent physical interpretation. The defining equation for G(z) 
Eq. (|2.8| ) implies that G(z) is an analytic function with two branch cuts between 
[—a + |,Q! + |] and [—a — |,a — |]. The matrix model saddle point equation 
identifies the values of G(z) along the two cuts via 

G(A + f±ie) =G(A-i=F«e); A G [—a, a]. (2.10) 

Thus the saddle point equation may be viewed as a 'gluing' condition for the two 
cuts on the z-plane. In addition G(z) = G(z*)* and G(z) = G(—z). The function 
G(z) is then uniquely determined by its asymptotic behaviour at z — » oo which is 
specified by the polynomial U(z) in Eq. ( |2.8|) . 

Given the density of eigenvalues p(X) in the matrix model one may use it to com- 
pute expectation values of various quantities in the matrix model. However, at the 
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Figure 1: The z-plane and its pre-image in the u-plane. The cuts in the z-plane are 
indicated by the dotted lines with Roman letters to dictate how they are mapped from the 
u-plane. Three representative closed contours Ca, Cb and Ca> are illustrated. 



outset we would like to point out these expectation values are distinct from vacuum 
expectation values of the corresponding observables in field theory. In particular, 
expectation values of the quantities Tr $ n in the matrix model are defined in the 
obvious way: 

/a 
p(X)X n d\. (2.11) 
■a 

Using the formulae above one can relate these expectation values to the terms in the 
large- z expansion of G(z). However these expectation values ((Tr$ n )) are not to be 
identified with the corresponding condensates (Tr$ n )in the field theory, 

(Tr$% FT ^((Tr$")). (2.12) 



2.3 Emergence of the r-torus 



The two-cut structure of the z-plane (from the analytic properties of G(z)) and the 
gluing conditions Eq. ( |2.10| ) uniquely specify a torus which we denote by E? where 
f is defined to be the complex structure parameter of this torus. The contour Ca 
in the z-plane enclosing the cut [—a + + ~] maps to the A-cycle of the torus 
while the contour Cb joining the two cuts maps to the -B-cycle of the torus. We also 
define Ca> as the contour enclosing the lower cut. Note that Ca + Ca> can be pulled 
away from the cuts to enclose the point at infinity. These contours in the z-plane 
are indicated in Figure 1. 

The torus Ef may be thought of as usual, as a quotient of the complex plane, say 
the w-plane (this auxiliary variable u will lead to the required elliptic parametriza- 
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tion), by a lattice with periods 2uo\ and 2^2 with Ct^A^i — t. The physical signifi- 
cance of this torus will become apparent upon implementing the procedure of [1-3] 
for computing the field theory superpotential. In fact f will be identified with r/N, 
signalling a low-energy SL(2, Z) duality symmetry in r/N where r is the coupling of 
the M = 4 theory. This duality, dubbed S-duality in [5] is distinct from the SL{2, Z) 
symmetry of the M = 4 theory although it is of course a consequence of the latter. 

The map z(u) from the it-plane to the z-plane is specified by the requirements 
that going around Ca (the A-cycle of the torus) returns z to its original value, while 
traversing the contour Cb (the S-cycle of the torus) causes z to jump by an amount 
i, which is the distance between the two cuts in the z-plane. Both these operations 
leave G unchanged which is therefore an elliptic function on the w-plane. Thus 

A - cycle : z[u + 2wi) = z{u); G(z(u + 2^)) = G(z(u)) , (2.13a) 
B- cycle: z(u + 2u 2 ) = z{u) + i; G{z{u + 2u 2 )) = G{z(u)). (2.13b) 



These requirements (along with the power law behaviour of G(z) at large z) 
specify z(u) uniquely and we find that 



z(u) = - 



7T 



C(u) U 



l^i(^l f ) 



2fi(^|f) 



(2.14) 



Here ((z) is the Weierstrass zeta function and 6\ is the Jacobi theta function (see [14] 
and the Appendix for details). Note that the map from the torus to the eigenvalue 
plane is completely independent of the choice of the field theory deformation V(§). 
Different deformations imply different asymptotic properties of the elliptic function 
G(z(u)) as z — > oo and hence lead to elliptic functions with different orders. 



2.3.1 An Example 

For the simplest deformation, i.e. the quadratic one corresponding to the N = V 
theory, the asymptotics determine G(z(u)) uniquely: 



V($) = $ 2 : G(z(u)) = ^ p{u) - 2^^ . (2.15) 



7T 2 
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G(z(u) ) will involve higher powers of the Weierstrass-p function for more complicated 
deformations. 

Let us see how in the simplest example which was first discussed and solved 
in [3,4], V(<&) = $ 2 , f gets identified with r/N. To see this we integrate the one- 
form G(z)dz over the A-cycle and 5-cycle to obtain the total number of eigenvalues S 
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and dJ-'o/dS respectively, the latter being the variation in the genus zero free energy 
of the matrix model upon transporting an eigenvalue to infinity, 



114 = 2mS = -i [ G(z(u))^du = i^P, (2.16a) 
J A du 12 GST 

„ dTn f ^, , ..dz(u) , f dEiif) 1 „ . . 

Then the effective superpotential W e s for the M = 1* theory in the confining vacuum 
is obtained by extremizing the following expression with respect to S, 

W eS = NTl B - tTI a . (2.17) 

Using our expressions for the respective integrals over the cycles it is easy to see that, 

and 

W eS = -^E 2 (t/N) . (2.19) 

Upon restoring the mass scales in the problem this result agrees precisely, up to 
an additive constant with the expressions of [6,7]. The fact that S in Eq.( [2.16a| ) 
automatically comes out as the r-derivative of the effective superpotential is in line 
with the identification of S with gluino condensate of the gauge theory. 



The identification of the complex structure f of the matrix model torus with 
t/N in the confining vacuum of the field theory in fact holds for all possible choices 
of V($) and indeed, seems to work for a completely different class of deformations 
namely the Leigh-Strassler deformations which will be discussed in detail below. The 
emergence of a torus with complex structure t/N, controlling the low-energy physics 
of the confining vacua of deformations of the M = 4 theory, has been well understood 
from the point of view of Seiberg-Witten theory in [5,12]. It is quite remarkable to 
see this structure emerge from the matrix model. 



2.4 QFT eigenvalues and general superpotentials 

Before embarking on a discussion of the eigenvalues below, we should emphasize a 
subtle point that will emerge in due course. It will turn out that the distribution of 
eigenvalues in the matrix model described by the density p(X) is not to be identified 
with the distribution of eigenvalues of $ in the field theory at large N. In fact, the 
main result of this section is that the Dijkgraaf-Vafa proposal is tantamount to iden- 
tifying the eigenvalues of the adjoint scalar in the confining vacuum of the field theory 
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with the function z(u) in Eq. ( |2.14j) (after a simple shift of its argument), distributed 
evenly over a cycle on the E? torus. In contrast the matrix model eigenvalues are 
not uniformly distributed along this cycle. 

The function z(u) that we have obtained above is a map from the torus E? to 
the plane of eigenvalues z. In particular, z(u) maps the A cycle of the torus to the 
cut enclosed by Ca, 

C A * z(2uj iX + u 2 ) = "2^^y + $ x ^ H §1- ( 2 - 20 ) 

This cut extends from —a + ~ to a + | and represents the possible values of the 
eigenvalues A of $ shifted by i/2 (the i/2 shift is deduced simply from Eq. ( |2.8|) ). In 
view of this we define the function 

w s 1 8'n(7Tx\f) r 1 1 i / \ 

AW = -2») ; iehU '- <2 ' 21) 

At first sight, this is simply a map from the A cycle of the torus to the interval 
[—a, a] in which the eigenvalues are distributed. However, we will soon find that 
X(x) has greater significance. Before we uncover this, we note that given the map 
X(x) we can ask how the eigenvalues of the matrix model are distributed along the 
A cycle of the torus. The answer is rather complicated being given by 

dn , , . ..d\ , s 

s - P(A(I % • (2 - 22) 

with — xq < x < Xq where X(xq) = a. 4 

In contrast, we will now show how the proposal of [4] automatically implies that 
eigenvalues of $ in the confining vacuum of the field theory at large N are uniformly 
distributed over the x-interval [— ~, |] and are given precisely by the function X(x). 
To see this we simply note that the function z"{u) = —p'(u) is elliptic on the torus 
and thus the following integral vanishes 



>C A +C' A 

Using the definition of G(z) in Eq. ( p.8|) this implies 



/ G{z(u))z"{u)du = 0. (2.23) 



S (uj(z+^)-uj{z-t))z"{u)du = i U(z)z"(u)du. (2.24) 

Jc A +C' A " Jca+C' a 

The periodicity of z(u) along the A-cycle allows us to reduce the left hand side of 
this expression to a contour integral around u = picking out the residue at u — 0. 
4 Notice that the A cycle covers the interval [—a, a] twice, so one must choose one of the branches. 
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Using z(2uix ± lu 2 ) = X(x) ± |, and Eq. ( |2.9| ) and the heat equation for the Jacobi 
theta functions (|A.9|) we find that 



1 d f l/2 

S = -W-^LI V{\{x))dx. (2.25) 

-1/2 



2m dr 



with A(x) given by Eq. ( |2.21| ) Since the proposal of [4] implies the identification 
f = t/N and 2niS = —dW^/dT so that 5 is the gluino condensate, we must have 

W eS = N V (A(x)) dx- X(x) = -- 3 ' . (2.26) 

7-1/2 2 6 3 {7TX\T) 

This proves that the eigenvalues of the adjoint scalar $ in the confining vacuum 
of the field theory at large N are given by the function X(x) and are uniformly 
distributed over the interval x G [— |, |] independent of the deformation V(3>). Note 
that, on the contrary the distribution of eigenvalues of the matrix model depends on 
the deformation in a complicated way. 

We should also emphasize that the result (|2.26| ) is also, according to Dijkgraaf 
and Vafa, true at finite N. At finite N the eigenvalues of $ in the confining vacuum 
have been computed from a completely different approach using the connection be- 
tween M = 1* theories and integrable systems in [8]. The eigenvalues found in [8] 
depend on a discrete index j and were given as 

Aj 2 9 3 (n/2-nj/N\r/Ny 3 A ' ' ' ' ' ' ^ 

(Notice that these, as explained above, are not to be identified with the eigenvalues 
Xj of the matrix model!) The Dijkgraaf- Vafa prediction from the matrix model above 
in Eq. ( |2.26| ) agrees precisely with the results of [8] upon taking the large-iV limit 



where the discrete label j in Eq.( p.27|) is replaced by the continuous variable x. In 
particular the eigenvalues are uniformly distributed in x. This demonstrates that 
using the recipe of [1-3] the matrix model indeed computes the superpotentials for 
the entire class of deformations considered above, provided N is large. Alternatively, 
for any given deformation it also computes all the condensates (Tr $ n ) via 

(Tr$ n ) = N X n {x)dx . (2.28) 

J-l/2 

In concert with all our previous warning (Tr<3> n ), the condensate of the field theory 
is not to be identified with the expectation values in the matrix model ( |2.11|) : 

((Tr$ n )) ^ (Tr$ n ) . (2.29) 
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As we have already mentioned, the eigenvalues and hence all the condensates in 
the confining vacuum are completely independent of the actual choice of deforma- 
tions. This is in complete accord with expectations from field theory. The defor- 
mations of the Af = 4 theory that we have discussed may be thought of as Af = 1 
perturbations of the Af = 2* theory. The confining vacua of the resulting theories 
arise from points on the Coulomb branch where the Seiberg-Witten curve for the 
Af = 2* theory degenerates maximally into a torus of complex structure r/N. The 
location of such a point in moduli space given by condensates (Tr$ n ), is completely 
independent of the choice of V($). This is in fact expected to be true in all massive 
vacua of the Af = 1 theory. 

Although we have seen a remarkable agreement between the matrix model recipe 
and the field theory results of [8] , strictly this agreement appears when the eigenvalues 
form a continuous distribution and can happen only in the large- N limit of the field 
theory. Indeed, the appearance of integrals over a continuous label rather discrete 
sums in Eq. ( p. 26 ) suggests that we are really doing a large- N calculation from the 



point of view of the field theory. So does the large- N matrix model, or the planar 
graphs of the matrix model really compute the holomorphic sector of the finite-iV 
field theory? The resolution of this puzzle will be the subject of Section 3. 



2.5 Some explicit examples 



We have computed the superpotentials from the matrix model by explicitly obtaining 
G(z(u)) for V(<&) = $ 2 and $ 4 and computing S, dJ^o/dS as the integrals of G(z)dz 
over the A and B cycles. These results match with a direct computation of the 
superpotential using the eigenvalue distribution following Eq. (|2.26|) . We quote the 
results for these cases below 

(i) \/($) = $ 2 : 

G(z(u)) = ^p(«) - *^ , (2.30a) 

n. = %); Z^-m. (2.30b) 
W cS = -Nh(r/N) = -^E 2 (t/N). (2.30c) 

The additive normalization of W e g can be fixed by a direct computation of Eq. ( [2.26] ) 
which gives 

r 1 / 2 N 
W cS = N \\x)dx = -—(E 2 (r/N) - l). (2.31) 

J-l/2 lZ 
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(ii) V($) = $ 4 



G(*(u)) = 4p(u) 2 + (vr 2 - 8^CK)) 

- + 6tt 2 C(wi) - 36 Wl <(wi) 2 ) , (2.32a) 

n , = m.. n B = f d -iP- Mr). (2.32b) 

W eff = -Nh(r/N) = ^-(10E 2 (t/N) 2 - E 4 (r/N) - 30E 2 (r/N)). (2.32c) 

The additive normalization of W e s can be fixed by a direct computation of Eq. fl2.2G| ) 
which gives 



r 1/2 n 

W cS = N A 4 (x) dx = —(10E 2 (t/N) 2 - E 4 (t/N) - 30E 2 (r/N) + 2l). (2.33) 
7-1/2 < Z(} 



2.6 Summary 

In summary, in this section we have verified that the Dijkgraaf-Vafa prescription 
for computing M = 1 superpotentials from matrix models works for a large class of 
deformations of the M = 4 theory. In particular it appears to provide a recipe for ex- 
tracting quantum field theory eigenvalue distributions from the matrix model which 
therefore permits an exact computation of all higher condensates and superpoten- 
tials for general deformations of the type discussed above. Our check of this proposal 
rests on the precise agreement between these field theory eigenvalues extracted from 
the matrix model, and those computed independently from a completely different 
approach in [8]. 

It must be emphasized that the field theory eigenvalue distributions are not 
identical to the matrix model distributions. The spectral densities of the two systems 
are different although there exists a recipe for extracting the QFT distributions from 
those of the matrix model. 



3. Large N versus finite TV 



As we have mentioned, one remarkable consequence of the Dijkgraaf-Vafa hypothesis 
is that the exact superpotential is captured by a large- N analysis. This presents a 
real puzzle: we have already remarked that the matrix model yields the distribution 
of eigenvalues of $ which agrees with the large- N limit of the distribution computed 
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from the integrable system approach in [8]; however, surely when we calculate the 
condensates (Tr $ n ) using the finite- iV distribution of [8] , we will find differences 
from the large- N result. 

Put another way, the large-iV matrix model (DV) result for the condensates 
( |2.28| ) are given by an integral along one of the cycles on the torus which can be 
written in terms of an integral over theta functions using Eqs. ( |2.21 ) and (|2.28j ): 



p^-srUS&tMi ix , (3.1; 



_ 1/2 V 2 3 (itx\t/N), 

while, in the same normalization, the condensates at finite- iV can be calculated from 
the exact eigenvalues that follow from the integrable system. In the confining vacuum 
these are given in Eq. ( |2.27| ). Hence, the condensate which follows from the integrable 
system analysis is 

(Tr*»y -VV '«ft!»2V (3 2) 

Clearly, (|3.1| ) and (|3.2|) differ at finite iV and so, one might think, that at finite iV 
the whole Dijkgraaf-Vafa hypothesis is called into question. 

However, it is well-known that the condensates of the M = 1* theory are plagued 
by mixing ambiguities [5,15]. There is no canonical definition of the operator $ n 
and different approaches can lead to expressions which differ by admixtures of <3> p , 
p < n (including the identity operator p = 0). What is important, however, is 
that the mixing coefficients should be vacuum independent. This means that the 
mixing coefficients should only depend on modular functions involving the torus 
whose complex structure is the bare coupling r rather than the effective coupling in 
the confining vacuum r/N. It simply would not make sense if the mixing coefficients 
between two approaches depended on the vacuum — in the present case the confining 
one. Having appreciated the possibility of mixing we can now see how the Dijkgraaf- 
Vafa hypothesis can be correct even though fl3.1|) and (|3.2|) are manifestly different 
at finite N. The resolution is that definitions of (Tr$ n )Dv and (Tr$ n ) int differ by 
vacuum independent mixings. 

Although we have no general proof of this, we can present very strong arguments 
for its verisimilitude. Firstly consider n = 2. In this case, the sum ( |3.2| ) gives 

(Tr $ 2 ) int = -?L(E 2 (t/N) + (N — 2)E 2 (r) - (N - l)£ 2 (r/2)) . (3.3) 

Notice that ( |3.3|) actually vanishes for N = 2 which follow from the fact that the 
two eigenvalues vanish in this case. Notice also that it differs from ( |2.31| ) by terms 
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which involve modular functions of r and r/2, rather than t/N and, therefore, one 
suspects that they differ by a vacuum independent constant which corresponds to 
mixing with the identity operator. 

At this point we should remind the reader that the class of theories we are 
considering have a large number of vacua with a mass gap, where the theory is realised 
in different phases. The canonical examples are the Higgs vacuum where the gauge 
group is broken completely by the Higgs mechanism, and the confining vacuum where 
classically the full gauge symmetry is unbroken although quantum mechanically there 
are no massless states due to strong coupling effects which cause the elementary 
degrees of freedom to confine and generate a mass gap. It is well-known that the 
Higgs and confining vacua are actually related by an 5-duality, r — > — 1/r, while 
more general modular transformations map to other massive vacua [6,7,17]. This 
provides us with a natural definition for the condensates. 

One way to disentangle the operator ambiguities is to note there exists a natural 
definition of the condensates, which we denote simply as (Tr$ n ), that have "good" 
modular properties. What we mean by this is that the value of the condensate in 
the Higgs vacuum is given by taking the value in the confining vacuum with the S- 
duality replacement r — > — 1/r and then by multiplying by the factor r~ n [5,6]. The 
values in other vacua are also related by other modular transformations. That such 
a basis of operators with good modular properties should exist follows ultimately 
from the fact that we are considering a deformation of the M = 4 theory, which is 
self-dual under modular transformations, by operators which can be chosen to have 
a definite weight. By assigning a compensating modular weight to the couplings 
the resulting deformed theory is then modular invariant. However, the vacua of 
the theory spontaneously break the modular symmetry. The modular symmetry 
then simply permute the vacua. We can also establish the existence of the modular 
covariant basis of operators directly from the integrable system: this tells us how 
(Tr<£> n ) is related to (Tr$ n ) int and therefore, using (|3.3| ) and its generalization to 
arbitrary deformation, how (Tr$ n ) is related to (Tr$ n )Dv- 

For the $ 2 operator, the modular covariant definition of the condensate is related 
to the matrix one by the vacuum independent — so involving modular functions of r 
rather than t/N — shift 



This then yields the modular covariant value of the condensate derived in [6] in the 
confining vacuum: 




(3.4) 



Confining: 




(3.5) 
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The expression for the condensate in the Higgs vacuum derived in [6] is 

N 2 

Higgs: (Tr $ 2 > = - — (NE 2 (tN) - E 2 (r)) . (3.6) 

As explained above, the expressions ( |3~5| ) and ( |3.6| ) are related by 

(Tr$ 2 ) = T~ 2 ((Tr$ 2 ) ) . (3.7) 

Higgs V Conf/r^-l/r 

In order to see this one uses the modular transformation property of the Eisenstein 
series (|ATJ). 

Now consider the $ 4 deformation. In this case the sum ( p.2|) gives 

(Tr $ 4 ) int ~{lOE 2 (r/Nf - E 4 (t/N)) 

+ i((2tf - 3)£ 2 (r/2) - (2N - 6)E 2 (r)) (Tr $ 2 ) int 

+ ^( 5N ( N ~ 1 )^(r/2) 2 - 10N(N - 2)E 2 (r) 2 
— (N — 1)(N - 3) 2 E 4 (r/2) + iV 2 (iV - 2)E 4 (r)) . 

In the large- N limit one can easily see that ( |3.8| ) goes over to ( f2.33j ). In addition it 
vanishes for iV = 2 as remarked above. Notice that the difference between (Tr $ 4 )dv 
and (Tr$ 4 ); nt is again by a vacuum independent, i.e. r/N- independent, admixture 
of Tr $ 2 and the identity operator. Once again there is a natural modular covariant 
definition of the condensate which can be established by investigating the integrable 
system and establishing a relation between (Tr$ 4 ) and (Tr$ 4 )i nt . Using (|3.8| ) one 
eventually finds 

(Tr$ 4 > = (Tr$ 4 ) DV + -^-(6-iV J E 2 (r))(Tr<I> 2 ) D v 

N 12 ( 3 - 9 ) 
+ — (N 2 E a (t) - 5N 2 E 2 (t) 2 - 35NE 2 (t) + 9) . 

The result ( [2.33j ), then gives the value of the condensate in the confining vacuum: 

Confining: (Tr $ 4 > = A (lOE 2 (r/N) 2 - E A {r/N) 

- 5N 2 E 2 (t) 2 + N 2 E 4 (r) - 5NE 2 (r)E 2 (r/N)) . 

The value in the Higgs vacuum follows by taking r — > — 1/r and multiplication by 
r~ 4 yielding 

Higgs: (Tr $ 4 > = ^L(lON 4 E 2 (rN) 2 - N 4 E 4 (rN) 

- 5N 2 E 2 (t) 2 + N 2 E a (t) - 5N 3 E 2 (t)E 2 (tN)) . 
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4. Perturbation of the Leigh- Strassler deformation 



We now turn to an application of the techniques of Dijkgraaf and Vafa to a very 
different class of deformations of the U(N) Af = 4 theory. We consider the M = 4 
theory perturbed both by a certain exactly marginal operator (the Leigh- Strassler 
deformation [18]) and certain mass terms. (The authors of [4] have pointed out 
that such deformations can accessed via known matrix models.) The Leigh- Strassler 
deformations of the A~ = 4 theory lead to a 3-dimensional fixed manifold of Af = 1 
theories, via special trilinear perturbations of the Af = 4 superpotential. Among 
these marginal deformations, of particular interest is the so-called "g-deformation" 
of the Af = 4 superpotential, 

W = Tr(i$[$ + ,$"] /3 ) = Tr(i$[$+$-e i/3/2 - ^-^ + e~ ip/2 ]) . (4.1) 



We consider a special perturbation of this theory via the superpotential 

W = Tr(i$[$ + ,$-] l3 + $ + $- + $ 2 ) (4.2) 

We choose this particular form for the mass perturbation as the corresponding large- 
N matrix model can be solved [19] around the trivial classical solution <£> = $ ± = 0. 
This trivial solution corresponds of course to the confining vacuum of the Af = 1 
gauge theory. 



As demonstrated in [19], expanding around the trivial vacuum the fields $ ± 
can be integrated out and the corresponding one-matrix model can be solved in the 
large- N limit by going to an eigenvalue basis and a change of variables. We simply 
quote this change of variables here, without further explanation. In particular the 
eigenvalues Aj of $ are redefined in terms of new variables 5i 

The eigenvalues experience a force that causes them to condense and expand out 
into a cut in the large- N limit. The complex z-plane in which the <Vs live turns out 
to be a cylinder with two cuts and gluing conditions on the two cuts implied by the 
saddle point equation for the function G(z). 



4.1 The t torus 



The key point for us is that large- N matrix model solution once again leads to the 
emergence of a torus. The z-plane being a two-cut cylinder with certain gluing 
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conditions, there is a natural map to a torus E? with 2 marked points. As before, we 
think of this torus as the complex w-plane modded out by lattice translations. The 
map z(u) and G(z(u)) can be determined precisely using in particular the asymptotic 
behaviour of G(z(u)). We refer the reader to [19] for the details of these asymptotic 
properties. Using these properties we are able to determine G(z(u)) and the map 
z(u) in the elliptic parametrization as 

mi , u w i cos(/?/2) 9 1 ([3/2) 2 r . oln , . ot , A AS 

and 

/N 1 x 0UO) 0i(7tu/2u;i - 5/4) , , 

expz n = -cot 09/2 ^Vttt — i — - — ^ttt- 4.5 

Note that G(z(u)) is elliptic while z{u) is periodic only along the A-cycle while 
S-cycle shifts lead to shifts in z(u), 

z(u + 2u{) = z(u); z{u + 2uj 2 ) = z{u) +i/3. (4.6) 

This encodes the gluing conditions implied by the saddle point equation satisfied by 
G(z). Note that the points u = ±/3co>i/27r are the two marked points on the f-torus 
and they map to the points z = =Foo. Note also that in the (3 — > limit z(u) and 
G(z(u)) reduce to Eqs.fl2.14|) and (|2.15|) after a simple rescaling of z(u) by (3. 



4.2 The effective superpotential 



We can now use the Dijkgraaf-Vafa prescription to compute the effective superpoten- 
tial for the relevant deformation of the Leigh Strassler theory. We have to compute 
the integrals of the one-form G(z)dz over the A and B cycles of the f-torus. The 
integrals are surprisingly simple and we find 

n A = -i [ G(z(u))z'(u)du = ^P-; IL B = -i [ G{z{u))z'{u)du = f^^--h(f) 

J A « r JB « r 

(4.7) 

where 

h(f) - COS{(3/2) 6limf) (4 8) 

n[T) 4sin 3 (/3/2) 9[(l3/2\fy ^ 

Thus we find the same structure as before for S and dT§ / dS which guarantees that 

dW cS /dS = t = t/N (4.9) 

and 

W - Nh - N COs(/3/2) 9 lM3lM (a in) 

WeS - ~ Nk - -iV 4simW2) W/2\r/N) ^ 
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up to an additive constant. A quick check confirms that as (3 — > this result 
reproduces the superpotential for J\f = 1* theory. 5 



4.3 QFT eigenvalues and condensates 



We recall that in the deformations considered in Section 2, the map from the A-cycle 
of the the torus E f to the 2-plane turned out to yield the field theory eigenvalues 
distributed uniformly on this cycle of E?. We can ask if a similar picture emerges 
for perturbations of the Leigh- Strassler theory. We first note that we need to take 
into account the field redefinition Eq. ( |4.3|) and so we need to consider the function 
— exp z{u) + l/(2sin(/?/2)) evaluated along the A-cycle of the torus, and we find that 
the relevant function (after a shift of z by i/3/2) is 

1 /ax 0UO) fl 3 (7ra-/3/4) 1 . -, i , 

^ = mkelJm + ^m xe[ ^ v (4 ' n) 



It is quite remarkable that for the quadratic deformation considered in this sec- 
tion we find that indeed, 

w M = nJ^x = -ivj-^jy ewWm + (4 ' 12) 

Although we don't have a general proof, the above leads us to believe that X(x) 
indeed yields the eigenvalues of the adjoint scalar in the confining vacuum of the field 
theory, distributed uniformly over the interval [— In other words, replacing x 
by a discrete index as in the previous section would precisely yield the eigenvalues 
of the adjoint scalar in the confining vacuum of the finite N field theory. We further 
conjecture that our conclusions in Section 3 on the relation between the field theory 
quantities at large- N and finite iV should also hold for the deformations of the Leigh- 
Strassler theories. Note also that our earlier comments on the distinction between 
field theory eigenvalue densities and matrix model eigenvalue densities apply in this 
case as well. 



A potential test involving S -duality: 

Although not completely understood in the context of this class of theories, 
one expects the Higgs and confining phases of these theories to be related by S- 

5 Actually there is a r-independent additive constant which diverges in the f3 — > limit which 
is a consequence of the linear shift of variables introduced in [19]. Strictly speaking, the resultant 
shift in the matrix model action 7V/(4sin 2 (/3/2)) needs to be reintroduced into W g which would 
then get rid of the offending additive term in the (3 — > limit. 
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duality i.e. r — ► — 1/r. A possible non-trivial test of our results for the superpo- 
tential and eigenvalue distributions in the confining vacuum is a comparison with 
classical eigenvalues in the Higgs vacuum after performing 5-duality on the confin- 
ing vacuum and taking the semi classical limit. It is not difficult to see that the 
classical eigenvalues in the Higgs vacuum are distributed on a circle with Xj ~ 
exp(-i(N + l-2j)f3/2)/cos(N/3/2); j = l,...N. After an ^-duality transformation, 
and redefining (3' = t/3, in the semiclassical r — > ioo limit our result Eq.( |4.11| ) also 
yields a uniform circular distribution of eigenvalues X(x) ~ exp(— i(3'Nx)/cos(f3N /2); 
x G [—1/2, 1/2]. This gives a hint of how ^-duality should act on the Leigh- Strassler 
deformation parameter. Indeed the Leigh-Strassler family of M = 1 CFTs are not 
self-dual under SX(2,Z) [18]. Our observations here could provide a starting point 
for further investigation into these issues. 



5. Conclusions 

In this article we have performed a test of the proposal of [1-3] for a large class of 
M = 1 preserving deformations of the M = 4 theories in their confining phase. We 
have shown that the Dijkgraaf-Vafa proposal translates into a rather simple recipe 
for extracting eigenvalues of the adjoint scalar $ in the confining vacuum from the 
corresponding matrix model, thereby yielding all the relevant condensates including 
the superpotential in that vacuum. We find that the eigenvalues extracted using this 
recipe match precisely with the results of [8] wherein the relevant eigenvalues were 
obtained from a completely different approach involving the connection of M = 1* 
theories with integrable models. Although at first sight the matrix model recipe yields 
large- N field theory results, we demonstrate through certain nontrivial examples that 
these large-A^ results are simply related to their finite N counterparts via vacuum 
independent field redefinitions or operator mixings. Thus the large- N matrix model 
indeed computes the holomorphic observables of the field theory at finite N in its 
confining vacuum. 

We also compute the effective superpotential in the confining vacuum of a mass 
deformation of the Leigh-Strassler CFT. We demonstrate that for this class of de- 
formations, the matrix model once again yield the field theory eigenvalues and con- 
densates via a simple recipe. A potential test of our results in the confining vacuum 
involves utilising S'-duality to obtain the eigenvalues in the Higgs vacuum in the 
classical limit. The latter can be computed independently and we find nontrivial 
agreement if we assume that the Leigh-Strassler deformation parameter transforms 
with modular weight one, under S-duality. 
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An interesting question concerns whether the matrix model approach can capture 
the full vacuum structure of the Af = 1* theory, including the Higgs vacuum. We 
answer this in the affirmative in a separate publication. 
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AppendixA: Some Properties of Elliptic Functions 

In this appendix we provide some useful — but far from complete — details of ellip- 
tic functions and their near cousins. For definitions and a more complete treatment, 
we refer the reader to one of the textbooks, for example [14]. An elliptic function 
f(u) is a function on the complex plane, periodic in two periods 2ui and 2uj 2 - We 
will define the lattice T = 2uj{L © 2072^ and define the basic period parallelogram as 

V = {u = 2/ilu 1 + 2vlu 2 , < n < 1, 0<z/<l}. (A.l) 

The complex structure of the torus defined by identifying the edges of V is 

t = cj 2 M (A.2) 

and we also define 

q = e inr . (A.3) 



A.l The Weierstrass function 



The archetypal elliptic function is the Weierstrass p(u) function. It is an even func- 
tion which is analytic throughout V, except at u — where it has a double pole: 

p{u) = — + V c k+1 u 2k , 

(A.4) 

92 93 3 \^ 1 ^ A 

C2= 20' C3= 28' Ck = { 2k + l)(k-3)^ Ck -> k ~ 4 - 
The Weierstrass function satisfies the fundamental identity 

(^PH) 2 = 4p(«) 3 - g 2 p(u) - g 3 (A.5) 
which defines the Weierstrass invariants #2,3 = #2,3 {wt) associated to the torus. 
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A. 2 The Weierstrass zeta function 



We are also interested in other functions which are only quasi-elliptic. First we have 
C(u). It is an odd function with the quasi-elliptic property: 



C(u + 2^) =CH + 2CM . 

Its derivative gives minus the Weierstrass function 

p{u) = -C(u) . 

It follows that has a simple pole at u — 0. Useful identities are 

u 2 CM - WiCM = y » 

C(wi + cj 2 ) = CM + CM , 

p(ui + uj 2 ) + pM + p(u 2 ) = . 



(A.6) 



(A.7) 



(A.8a) 

(A.8b) 
(A.8c) 



A. 3 The Jacobi theta functions 



We are also interested in the Theta functions 9i(x\r), or 9i(x,q), i = 1,2,3,4. They 
are also quasi-elliptic functions on V when x = ttu/2uji. Each of them satisfies the 
heat equation 

(A.9) 



d 2 9AxW) .dOAxW) 
7ri — + 4 V = . 



dx 2 Or 

They are related to the previous functions; for instance, two relations that we need 
are 



C(u) U 



7T 9' 1 (x\r) 



2ui 9\{x\r) 



in 7T 9' 3 (x\t) 



(A.10) 



where the derivative is with respect to x. 



A. 4 The Eisenstein series 



Finally, we introduce the Eisenstein series and their relation to elliptic functions. 
There are a number of ways to define these series: see [16] 6 ; for instance 

EJt) = - V — -7 , (A.ll) 

m.neZ v ' 
(m f nj — 1 



6 Beware, the definition of q in that reference is the square of ours. 
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where (to, n) denotes the greatest common divisor. Each series has a g-expansion; 
for example 

n=l 

E 4 (r) = l + 240^a 3 (n)g 2 ", (Al2) 

n=l 

E 6 (r) = l- 504j2^(n)q 2n , 

n=l 

where <Jj{n) is a sum over each positive integral divisor of n raised to the j th power. 
Useful identities are 

E 2 — -q i E i ~ E *) , 

Ei= 2 f(E 4 E 2 -E 6 ), (A-13) 

E' % = m{E & E 2 -El) , 

where the derivatives are defined with respect to r. Under a ^-duality transformation 
we have 

E 2 (-1/t) = r 2 E 2 (r) + — , E k (-l/r) = r k E k (r) k > 4 . (A.14) 

The relations of the Eisenstein series to the elliptic functions defined previously 
appears through the identities 

7T 2 4-7T 4 8-7T 8 

CMfi = ^^(r) , ^ 2 (2 Wl ) 4 = ^-^4(r) , s 3 (2^) 6 = ^^ 6 (r) . (A.15) 



AppendixB: Stationary points of the Hamiltonians 

In this Appendix we show that the configuration of the ECM system considered in the 
text is an equilibrium point with respect to each of the Hamiltonians H p = TyL p . In 
each case Hamilton's equations read dPj/dt = —dH p /dXj and dXj/dt = dH p /dPj. 
Thus we need to show that the right hand side of both these equations vanishes 
in the configuration of interest where Xj = (Xj) = 2mj/N for j = 1,2... iV and 
Pj = (Pj) = 0. (Note that here for simplicity we present the requisite proof for the 
equilibrium positions that correspond to the Higgs vacuum of the gauge theory. The 
proof for the confining vacuum can be obtained along the same lines by using the Lax 
matrix appropriate for that vacuum and equilibrium positions Xj = (Xj) = 2nirj/N 
for j = 1, 2 ... N and Pj = (Pj) = 0.) 
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Importantly we must also remember to impose the constraint appropriate to the 
SU (N) theory: H\ = X^Li Pj = 0- This is most easily accomplished by introducing a 
Lagrange multiplier ji and looking for stationary points of the modified Hamiltonians 
H p = H p — fJ,H 1 . Thus we must show that the equations, 



dHp 



= 



dHp 
dPi 



= V 



(B.l) 



Can be solved in the equilibrium configuration. We will first demonstrate that, 



dHr, 



= 



(B.2) 



The Hamiltonians of the integrable system can be expressed with the Lax matrix as 



H p = Tr \U 



(B.3) 



so we would like to prove that for the Calogero-Moser system 



dHp 



Tr 



.i 9L 
dX~ 



(B.4) 



The CM system is described with the Lax matrix 

L jk = PjSjk + i(l - Sj^xiXj - X k ) (B.5) 

where the function x(u) is given explicitly in [8]. For the present purposes we only 
need two properties of x{u). Firstly it has period 2ni: x{u + 2ni) = x(u). Secondly 
it is an odd function: x(—u) = —x(u). As above, the equilibrium positions of the 
integrable system are 

{Xi) = 2i^_ tI (N-l) j = 1 N (B6) 

and (Pj) = The key point is that the Lax matrix becomes a circulant when eval- 
uated in this configuration. It is then easy to show that the unitary transformation 
which diagonalises the Lax matrix is, 

L = U ] KU 

A jk = Xj8 jk (B.7) 



with 



v 
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The trace then can be written 

_! dL 



Tr 



L p ~ 



Tr 



A p ~ U 



(Xj),{Pj) 



dL 
dX~ 



(B.9) 



Substituting the expressions for the Lax matrix and for U the explicit form of 
the trace is 



Tr 



dL 



9 N N 9tt 

NY,*? 1 E ^{X 3 - Xl )^-{k-l){ 3 -l) 



k=l 



(B.10) 

Using that x'(u) is an even function and is periodic in u\ = 2iri we find the terms in 
the summation pairwise cancel. Thus the summation gives zero and we conclude 



dH n 



OX, 







(B.11) 



Similarly its easy to show that the second equation in ( B.l ) is solved by choosing 
/i = Hp-i. This completes the proof. 
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